Conformational Changes and Unfolding of ß-Amyloid Substrates in the Active Site of γ-Secretase.

Alzheimer's disease (AD) is the leading cause of dementia and is characterized by a presence of amyloid plaques, composed mostly of the amyloid-ß (Aß) peptides, in the brains of AD patients. The peptides are generated from the amyloid precursor protein (APP), which undergoes a sequence of cleavages, referred as trimming, performed by γ-secretase. Here, we investigated conformational changes in a series of ß-amyloid substrates (from less and more amyloidogenic pathways) in the active site of presenilin-1, the catalytic subunit of γ-secretase. The substrates are trimmed every three residues, finally leading to Aß40 and Aß42, which are the major components of amyloid plaques. To study conformational changes, we employed all-atom molecular dynamics simulations, while for unfolding, we used steered molecular dynamics simulations in an implicit membrane-water environment to accelerate changes. We have found substantial differences in the flexibility of extended C-terminal parts between more and less amyloidogenic pathway substrates. We also propose that the positively charged residues of presenilin-1 may facilitate the stretching and unfolding of substrates. The calculated forces and work/energy of pulling were exceptionally high for Aß40, indicating why trimming of this substrate is so infrequent.

Investigation of Peptides for Molecular Recognition of C-Reactive Protein-Theoretical and Experimental Studies.

We investigate the interactions between C-reactive protein (CRP) and new CRP-binding peptide materials using experimental (biological and physicochemical) methods with the support of theoretical simulations (computational modeling analysis). Three specific CRP-binding peptides (P2, P3, and P9) derived from an M13 bacteriophage have been identified using phage-display technology. The binding efficiency of the peptides exposed on phages toward the CRP protein was demonstrated via biological methods. Fibers of the selected phages/peptides interact differently due to different compositions of amino acid sequences on the exposed peptides, which was confirmed by transmission electron microscopy. Numerical and experimental studies consistently showed that the P3 peptide is the best CRP binder. A combination of theoretical and experimental methods demonstrates that identifying the best binder can be performed simply, cheaply, and fast. Such an approach has not been reported previously for peptide screening and demonstrates a new trend in science where calculations can replace or support laborious experimental techniques. Finally, the best CRP binder─the P3 peptide─was used for CRP recognition on silicate-modified indium tin oxide-coated glass electrodes. The obtained electrodes exhibit a wide range of operation (1.0-100 µg mL-1) with a detection limit (LOD = 3σ/S) of 0.34 µg mL-1. Moreover, the dissociation constant Kd of 4.2 ± 0.144 µg mL-1 (35 ± 1.2 nM) was evaluated from the change in the current. The selectivity of the obtained electrode was demonstrated in the presence of three interfering proteins. These results prove that the presented P3 peptide is a potential candidate as a receptor for CRP, which can replace specific antibodies.

GS-SMD server for steered molecular dynamics of peptide substrates in the active site of the γ-secretase complex.

Despite recent advances in research, the mechanism of Alzheimer's disease is not fully understood yet. Understanding the process of cleavage and then trimming of peptide substrates, can help selectively block γ-secretase (GS) to stop overproduction of the amyloidogenic products. Our GS-SMD server (https://gs-smd.biomodellab.eu/) allows cleaving and unfolding of all currently known GS substrates (more than 170 peptide substrates). The substrate structure is obtained by threading of the substrate sequence into the known structure of GS complex. The simulations are performed in an implicit water-membrane environment so they are performed rather quickly, 2-6 h per job, depending on the mode of calculations (part of GS complex or the whole structure). It is also possible to introduce mutations to the substrate and GS and pull any part of the substrate in any direction using the steered molecular dynamics (SMD) simulations with constant velocity. The obtained trajectories are visualized and analyzed in the interactive way. One can also compare multiple simulations using the interaction frequency analysis. GS-SMD server can be useful for revealing mechanisms of substrate unfolding and role of mutations in this process.

Specificities of Protein Homology Modeling for Allosteric Drug Design.

The allosteric binding sites are usually located in the flexible areas of proteins, which are hardly visible in the crystal structures. However, there are notable exceptions like allosteric sites in receptors in class B and C of GPCRs, which are located within a well-defined bundle of transmembrane helices. Class B and C evolved from class A and even after swapping of orthosteric and allosteric sites the central binding site persisted and it can be used for easy design of allosteric drugs. However, studying the ligand binding to the allosteric sites in the most populated class A of GPCRs is still a challenge, since they are located mostly in unresolved parts of the receptor's structure, and especially N-terminus. This chapter provides an example of cannabinoid CB1 receptor N-terminal homology modeling, ligand-guided modeling of the allosteric site in GABA receptor, as well as C-linker modeling in the potassium ion channels where the allosteric phospholipid ligand PIP2 is bound.

The Role of Cholesterol in Amyloidogenic Substrate Binding to the γ-Secretase Complex.

Alzheimer's disease is the most common progressive neurodegenerative disorder and is characterized by the presence of amyloid ß (Aß) plaques in the brain. The γ-secretase complex, which produces Aß, is an intramembrane-cleaving protease consisting of four membrane proteins. In this paper we investigated the amyloidogenic fragments of amyloid precursor protein (substrates Aß43 and Aß45, leading to less amyloidogenic Aß40 and more amyloidogenic Aß42, respectively) docked to the binding site of presenilin, the catalytic subunit of γ-secretase. In total, we performed 9 µs of all-atom molecular dynamics simulations of the whole γ-secretase complex with both substrates in low (10%) and high (50%) concentrations of cholesterol in the membrane. We found that, at the high cholesterol level, the Aß45 helix was statistically more flexible in the binding site of presenilin than Aß43. An increase in the cholesterol concentration was also correlated with a higher flexibility of the Aß45 helix, which suggests incompatibility between Aß45 and the binding site of presenilin potentiated by a high cholesterol level. However, at the C-terminal part of Aß45, the active site of presenilin was more compact in the case of a high cholesterol level, which could promote processing of this substrate. We also performed detailed mapping of the cholesterol binding sites at low and high cholesterol concentrations, which were independent of the typical cholesterol binding motifs.

Allosteric Modulation of the CB1 Cannabinoid Receptor by Cannabidiol-A Molecular Modeling Study of the N-Terminal Domain and the Allosteric-Orthosteric Coupling.

The CB1 cannabinoid receptor (CB1R) contains one of the longest N termini among class A G protein-coupled receptors. Mutagenesis studies suggest that the allosteric binding site of cannabidiol (CBD) involves residues from the N terminal domain. In order to study the allosteric binding of CBD to CB1R we modeled the whole N-terminus of this receptor using the replica exchange molecular dynamics with solute tempering (REST2) approach. Then, the obtained structures of CB1R with the N terminus were used for ligand docking. A natural cannabinoid receptor agonist, Δ9-THC, was docked to the orthosteric site and a negative allosteric modulator, CBD, to the allosteric site positioned between extracellular ends of helices TM1 and TM2. The molecular dynamics simulations were then performed for CB1R with ligands: (i) CBD together with THC, and (ii) THC-only. Analyses of the differences in the residue-residue interaction patterns between those two cases allowed us to elucidate the allosteric network responsible for the modulation of the CB1R by CBD. In addition, we identified the changes in the orthosteric binding mode of Δ9-THC, as well as the changes in its binding energy, caused by the CBD allosteric binding. We have also found that the presence of a complete N-terminal domain is essential for a stable binding of CBD in the allosteric site of CB1R as well as for the allosteric-orthosteric coupling mechanism.

The Hydrophobic Ligands Entry and Exit from the GPCR Binding Site-SMD and SuMD Simulations.

Most G protein-coupled receptors that bind the hydrophobic ligands (lipid receptors and steroid receptors) belong to the most populated class A (rhodopsin-like) of these receptors. Typical examples of lipid receptors are: rhodopsin, cannabinoid (CB), sphingosine-1-phosphate (S1P) and lysophosphatidic (LPA) receptors. The hydrophobic ligands access the receptor binding site from the lipid bilayer not only because of their low solubility in water but also because of a large N-terminal domain plug preventing access to the orthosteric binding site from the extracellular milieu. In order to identify the most probable ligand exit pathway from lipid receptors CB1, S1P1 and LPA1 orthosteric binding sites we performed at least three repeats of steered molecular dynamics simulations in which ligands were pulled in various directions. For specific ligands being agonists, the supervised molecular dynamics approach was used to simulate the ligand entry events to the inactive receptor structures. For all investigated receptors the ligand entry/exit pathway goes through the gate between transmembrane helices TM1 and TM7, however, in some cases it combined with a direction toward water milieu.